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Abstract. We present a database of 908 spectra of 709 stars 
obtained with the ELODIE spectrograph at Observatoire de 
Haute-Provence. 52 orders of the echelle spectra have been 
carefully fitted together to provide continuous, high-resolution 
spectra in the wavelength range AA = 410 — 680 nm. The 
archive provides a large coverage of the space of atmospheric 
parameters : T cff from 3700 K to 13600 K, log g from 0.03 to 
5.86 and [Fe/H] from -2.8 to +0.7. At the nominal resolution, 
R=42000, the mean signal-to-noise ratio is 150 per pixel. The 
spectra given at this resolution are normalized to their pseudo- 
continuum and are intended to serve for abundance studies, 
spectral classification and tests of stellar atmosphere models. 
A lower resolution version of the archive, at R= 10000, is cal- 
ibrated in physical flux with a broad-band photometric preci- 
sion of 2.5% and narrow-band precision of 0.5%. It is scoped 
for stellar population synthesis of galaxies and clusters, and for 
kinematical investigations of stellar systems. 

The archive is distributed in FITS format through the HY- 
PERCAT and CDS databases. 

Key words: atlases - stars: abundances - stars: atmospheres - 
stars: fundamental parameters - Galaxies: stellar content 



1. Introduction 

Spectral libraries covering the HR diagram with medium to 
high resolution and large spectral range are essential tools in 
astronomy. We have identified two areas where there will be a 
special need of such libraries in the coming years : the auto- 
mated parameterization of stellar spectra and the spectral syn- 
thesis of stellar populations of galaxies. New multi-object spec- 
trographs on large telescopes, like Giraffe on the VLT for ex- 
ample, will soon make it possible to observe hundreds of ob- 
jects during the same exposure with an unprecedented resolu- 
tion and spectral coverage. The challenge will be to extract the 
maximal information in a reasonable time. 

For stars, medium to high resolution spectra are the ba- 
sis for determining radial velocities, atmospheric parameters 
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(T g , log g , [Fe/H]) and eventually projected rotational veloc- 
ities, vsin i . A new method, called TGMET, is emerging to 
make the parameterization of stellar spectra in terms of tem- 
perature, gravity and metallicity fast and automatic (Katz et al 
|l998 ). It relies on the direct comparison of a target spectrum to 
a library of referenc e spec tra with known atmospheric param- 
eters (Soubiran et al 1998 ). TGMET was originally dedicated 
to the measurement of F, G, K stars observed with ELODIE. 
The present paper, which continues the TGMET project, par- 
ticipates to release these limitations. The two major points ad- 
dressed are (1) the extension of the reference library to all 
spectroscopic types and luminosity classes and its densifica- 
tion, and (2) to resample the spectra in wavelength (i.e. remove 
the instrumental signature in ELODIE spectra and connect the 
echelle orders together) in order to enable comparison of spec- 
tra from any origin. 

The TGMET project appeared to be also very useful for ex- 
tragalactic studies where medium resolution libraries with large 
coverage in stellar parameters are needed. Stellar libraries, cal- 
ibrated in flux, are used to model composite spectra through 
stellar population s ynthes is (eg. the PEGASE program, Fioc & 
Rocca-Volmerange 1997 ). An important issue is to disentangle 
the effects of age and metallicity and it has been demonstrated 
that th e spectral resolution is a k ey fact or (Worthey & Ottviani 



1997, hereafter W097, Vazdekis 1999). To take profit of a res 



olution higher than about R=2000 both the internal kinematics 
and the characteristics of the stellar populations must be fitted 
simultaneously. Making available an extended library of stellar 
spectra is the first step in this direction. 

Current libraries are based on low resolution spectroscopy 
(eg Jacoby et al. |1984[ Serote Roos et al |1996[ Pickles |1998| ) 
or they are restricted to a limited area of the HR diagram (eg 
Montes et al 2000 and Soubiran et al |l998| for F,G,K stars). 
Here we present an homogeneous dataset of spectra which fills 
the deficiency of such libraries. The sample of stars and associ- 
ated data are presented in Sect. 2. Sect. 3 & 4 describe the two 
phases of the data processing. In Sect. 5 we evaluate the qual- 
ity of the archive. The access to the data through HYPERCAlj] 



http://www-obs.univ-lyonl.fr/hypercat/! 1/spectrophotometry.html 
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or at the Centre de Donnees astronomiques de Strasbourg 
described in Sect. 6. 



2. Stellar content and data 

2.1. Physical parameters of the stars 

The present sample was assembled by me rging the stellar li- 
brary of F,G,K stars (Soubiran et al 1998 ) serving as a basis 
for the TGMET program with additional spectra taken from 
other ob serving programs and from the near-to-line ELODIE 
archive ( http://www.obs-hp.fi ). The spectra of the database cor- 
respond to stars which have been selected because they have 
published values of T c ff or (T e g , log g , [Fe/H]) , or reliable es- 
timations of the absolute magnitude My , deduced from Hip- 
parcos parallaxes. The main source of atmospheric parameters 
is the catalogue of [Fe/H] determinations, 1996 and 2001 ver- 
sions (Cayrel de Strobel et al 1997 and 2001) but two other 
sources have been also used : Carney et al ( 1994 ) and Thevenin 
( 1998 ). Effective temperatures were also taken in the lists pub- 
lished by Blackw ell & L ynas-Gray ( |1998[ ), Alonso et al ( |1996aj ) 
and Alonso et al (1999a), or calculated from the colour indices 



V-K or b-y following the relations T e g =f(colour, [Fe/H] ) es- 
tablished by Alonso et al 1996b| and 1999b for dwarfs and gi- 
ants respectively. Each star had (T c ff , log g , [Fe/H]) estimated 
by averaging the determinations found in the literature, giv- 
ing a half weight to old references and Stromgren photometry 
and a double weight to T c ff determined by Blackwell & Lynas- 
Gray ( 1998 ). As in the TGMET library, the parameters were 
then labelled according to their reliability which was estimated 
from the number of determinations and their standard deviation 
around the mean. The reliability scale runs from to 4, with the 
highest reliability 4 corresponding to uncertainties lower than 
80 K in T c ff , and 0.06 dex in [Fe/H] . These uncertainties can 
reach 1 15 K and 0.09 dex for reliability 3, and 150 K and 0. 1 1 
dex for the reliability 2. Reliability 1 is attributed for param- 
eters based on old determinations or photometry, or present- 
ing large discrepancies between references. Reliability means 
that no reference on atmospheric parameters was found in the 
literature. The database is made of 908 spectra corresponding 
to 709 different stars which are shown in the plane (log T ff , 
log g ) in Fig 1 and in the plane (log T G ff , [Fe/H] ) in Fig. 2. 
T c ff ranges from 3700 K to 13600 K, log granges from 0.03 to 
5.86 and [Fe/H] ranges from -2.8 to +0.7. In these two plots, 
new estimations of (T c g , log g, [Fe/H]) were used instead of 
those from the literature which are incomplete. These new es- 
timations were obtained with the current version of TGMET 
which is still under development (a paper presenting this new 
version is in preparation and a description is given in the elec- 
tronic version of the archive). Atmospheric parameters from 
the literature together with new estimations and measured line 
indices from the spectra presented here are given in Table 1, 
only available in electronic form (see sect. 6). 

More information on each star is available in the header of 
the FITS spectra. It includes the absolute magnitude My from 
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Fig. 1. T c ff vs. log g for the 709 stars of the archive. For the sake 
of homogeneity, the atmospheric parameters presented here are 
not from the literature which is incomplete, but were estimated 
with the current version of TGMET (Sect 5.3.3). 



Hipparcos parallaxe when available. A scale of reliability was 
also established from the precisi on of the parallaxe and V mag- 
nitude from Tycho-2 (H0g et al 2000 ), 4 corresponding to par- 
allaxes with a precision better than 10%. Also available in the 
FITS headers are the spectral types, (V, B-V) taken from the 
Tycho-2 catalogue (converted from the Tycho-2 Bt and Vt) 
and radial velocities usually measured at the telescope. The on- 
line cross-correlation technique was used to measure the radial 
velocity of strong-lined spectra, corresponding to moderate ef- 
fective temperatures up to 6500 K but for hotter stars, the ra- 
dial velocity was estima ted us ing a least-squares deconvolution 
technique (Donati et al. 1997 ). 

The database includes several stars with line profiles broad- 
ened by rotation, macroturbulence or binarity and a few stars 
with spectral peculiarities (Ap, Am, emission line stars...). 
These pecularities are indicated in the headers. 



2.2. ELODIE spectra 

ELODIE provides a spectral range of 390-680 nm recorded in 
a single exposure as 67 orders on a IK CCD at a mean re- 
solving power of 42000. Optimal extraction and wavelength 
calibration are automatically performed by the on-line reduc- 
tion software TACOS. Complete information on ELODIE can 



be found in Baranne et al (1996) and in the TACOS user man 



ual by Queloz (1990- The global efficiency of the instrument 



http://cdsweb.u-strasbg.fr/cats/III.htx 



drops in the blue by a factor 5 at 440 nm for a solar-type star 
and by a factor 50 at 390 nm. Hence we decided to limit the 
spectral range to AA = 410 — 680 nm. 
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Fig. 2. Distribution of the 709 stars of the archive in the plane 
T e ff vs [Fe/H] estimated with the current version of TGMET 
(Sect 5.3.3). 



The typical signal-to-noise (S/N) ratio per pixel at 550 nm 
is 150 (it is higher than 80 for more than 90% of the spectra). 
Several spectra with a lower S/N, down to 35, correspond gen- 
erally to faint deficient stars which have been included in the 
database for a better sampling of the parameter space. 

For the present purpose, the 908 spectra of the database 
were processed to provide ID spectra rebinned in wavelength 
and calibrated in absolute fluxes. Since the main goal of the 
archive was to construct a library of spectra representative of 
the various stellar types and an interpolated grid of spectra cov- 
ering the parameters space, all the spectra were reduced to the 
rest-frame. The archive was prepared from the flat-fielded and 
order-extracted TACOS spectra in two main steps. (1) The or- 
ders were de-blazed, the spikes due to cosmic rays and telluric 
lines were masked, all the orders were connected in one spec- 
trum, and the pseudo-continuum normalization was computed. 
(2) The spectra were calibrated into "physical" flux, ie. above 
the atmosphere. 



3. Strengthening and normalization to the 
pseudo-continuum 

This phase of the data processing takes the TACOS spectra 
to deliver order-connected spectra with a constant wavelength 
step. Both the instrumental and continuum fluxes are evaluated. 
The operations follow these steps: 

- Strengthen the spectra by dividing them by a de-blazing 
function 



- Compute the pseudo-continuum normalization and analyze 
the "shape" of each order to determine a correction to the 
diffuse-light subtraction. 

- Mask spikes and telluric lines. 

- Convolve and rebin with two desired resolutions (R= 10000 
and R=42000) to ID spectra. 

3.1. De-blazing and pseudo-continuum normalization 

The instrument, hard- and soft-ware, was mostly designed as 
a "velocimeter" and hence the spectrophotometric quality was 
not optimized. The standard flat-field correction and or der ex - 
traction was applied as described in Baranne et al. (1996). 



This procedure alone would not permit to connect the spectra 
since the blaze is not accurately enough corrected, each order 
presents a residual "curvature" of about 2%. This effect is cer- 
tainly due to diffuse light within the spectrograph which is im- 
perfectly corrected by the TACOS program, this is a common 
difficulty with echelle spectrographs. 

To solve this problem we found no other possibility but to 
analyze this residual curvature on the actual spectra and correct 
it. Obviously, it is a risky approach since the spectrum itself 
is used to determine the correction. To minimize possible bi- 
ases, the residual curvature is modeled with only a few param- 
eters : each order is represented by a second degree polynomial 
whose coefficients are not independent from order to order, but 
are themselves second degree polynomials of the order num- 
ber. The total effect of the diffuse light is modeled by only 3 
coefficients. 

Practically, we divide the TACOS spectrum by a de-blazing 
function determined from a tungsten-lamp spectrum (internal 
flat-field) and corrected using 3 metal-deficient stars (because 
the blaze function appears to have slight systematic differences 
between internal flat-field and star). 

To analyze the "residual" curvature of each order, we nor- 
malized the spectrum to its pseudo-continuum. The pseudo- 
continuum has been defined using one point in each order and 
then interpolated over the whole range of wavelengths with a 
3 rd degree spline. The pseudo-continuum level is determined 
by fitting a second degree polynomial to each order. The spec- 
trum is weighted with a mask intended to avoid as much as 
possible the stellar lines, the fit is iterated after k — a clipping 
and the remaining bias due to the weak stellar lines is statis- 
tically corrected. The latter correction is based of the analysis 
of the skewness of the distribution of residuals. The weighting 
mask is prepared from a previous iteration of the reduction of 
the whole archive: for each wavelength point it gives the prob- 
ability to have a stellar line (we started with no weight). 

Though the visual inspection of the strengthened spectra 
appear satisfactory, the detection of possible biases have been 
a major concern and is presented in Sect. 5. 

After the correction for diffuse light was applied, the 
pseudo-continuum was re-computed. These normalized spec- 
tra are scoped for abundance studies, and determination of the 
atmospheric parameters with TGMET. 
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All the information permitting calculation of the detected 
flux (in electrons) is kept in the form of a table extension at- 
tached to the FITS image. The S/N of each pixel is also kept in 
a separate extension. The detailed description of the FITS files 
is given with the electronic version of the tables and spectra. 

3.2. Cosmic rays and telluric lines 

The spikes due cosmics and defective pixels and the telluric 
lines have been m asked following the procedure developed by 
Katz et al. (|l998|). 



3.3. Resampling to ID spectra 

The orders are normalized to their pseudo-continuum and re- 
sampled to a ID spectrum with a wavelength step of 0.005 nm 
using spline interpolation (the observations are over-sampled 
by about 20%). 

To produce a version of the archive at the resolution 10000 
(i.e. 30 km.s _1 at A = 550 nm) the spectra are convolved by a 
Gaussian function of FWHM = 0.054 nm, and rebinned with a 
step of 0.02 nm. 

The high resolution version has a constant resolution of 
R=42000, ie. the resolution in nm changes with the wavelength. 
At variance, with the low resolution archive, the wavelength 
resolution, in nm, is constant throughout the spectra. 

Note also that the convolution to produce the low resolution 
archive is done on the continuum-normalized spectra, which is 
not formally equivalent to convolving the flux calibrated spec- 
tra. However the induced difference is absolutely insignificant. 

4. Flux calibration 

The goal of the flux calibration is to determine the spectral en- 
ergy distribution above the atmosphere (the physical flux). The 



standard procedure to achieve the objective (Bessell 1999 ) con 



sists in determining the instrumental response and correct for 
the atmospheric absorption by comparing the spectra to tem- 
plates of known spectral energy distribution observed with the 
same setup between the target stars. 

Since the primary scope of the observation was not spec- 
trophotometry, no particular care was taken in interleaving tem- 
plate stars and some observations were done through heavy at- 
mospheric absorption. Therefore, the standard procedure can- 
not be straightforwardly applied. Before entering into the de- 
tails of the procedure, we present hereafter its genesis. 

The relation between the instrumental flux deduced in the 
previous section and the physical flux (the flux calibration re- 
lation) primarily reflects the effects of (1) the instrumental re- 
sponse (ie. the combined effect of the whole optics and CCD 
spectral sensitivity), (2) the absorption by a clear atmosphere 
and (3) an additional extinction due to an atmospheric haze. 

We started from the assumptions that (1) the instrumental 
response is stable or at least stable during each observing run, 
(2) the absorption by the clear atmosphere can be parameter- 
ized by the airmass and (3) that the effect of the haze can be 



parameterized by a color excess defined as the difference be- 
tween the measured B-V color and the Tycho-2 B-V color. 

A fit of these three functions of the wavelength using spec- 
tra of stars with known spectral energy distribution (the de- 
tails are presented below) resulted in a photometric precision of 
about 5% when we considered the instrumental response stable 
over all observation runs. 

The shape of the three functions entering the calibration re- 
lation, and the analysis of the residuals lead to the suspicion 
that the correction of the atmospheric refraction provided by 
the spectrograph was not fully effective. This was diagnosed 
by ( 1 ) a stronger than expected wavelength dependence of the 
apparent atmospheric and haze transmissions and (2) a residual 
random "curvature" of the calibrated spectra that increases with 
airmass. Our interpretation is that the "clear atmosphere" func- 
tion corrects at the first order the effect of the refraction (the 
star is centered on the fiber aperture by the auto-guider on the 
blue end of the atmospheric spectrum). The uncertainty on the 
centering of the star (miss-guiding) results in a "color excess" 
absorbed by the haze function, but depending on the value of 
the miss-centering and of the seeing a residual curvature is left. 

It was unfortunately not possible to parameterize this effect 
since it depends on a random and unknown miss-centering and 
on the seeing which is not always recorded in the observing 
logs. An alternative would have been to use an additional color 
excess to constrain the curvature, but no such independent mea- 
surement is available for the stars of our archive. 

Adding a quadratic airmass term, another function of color 
excess and allowing for variation of the instrumental response 
between the different observing runs led to a photometric qual- 
ity of about 3% which was then improved to 2% after correc- 
tions based on internal comparisons. We did not find hints for 
other effect on the photometric quality, such as the effect of the 
stress on the optical fiber or positioning of the color correction 
filter that may change the instrumental response. 

We will now describe the different steps of the calibrations: 
(1) primary calibration based on the comparison with low spec- 
tral resolution flux templates (2) secondary calibration based 
on internal comparison and (3) correction of the instrumental 
response using template spectra with a better spectral resolu- 
tion than for the primary calibration. 

4.1. Primary calibration 

As already explained, the observations were not originally 
scoped for photometry and no flux templates were explicitly 
observed. However, since the archive contains bright stars, 
many of them are in common with other libraries or datasets. 
The largest intersection is with the collection of low resolu- 



tion spectra presented in Burnashev (1985, CDS III/126). The 
spectral resolution is about R=100 (FWHM = 5-6 nm), and 
the observations are from different sources which were homog- 
enized to a common scale. The observations were done with 
photo-electric spectro-photometers mostly a the Crimean Ob- 
servatory, and are not corrected from the interstellar extinction. 
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Some stars in our archive are in common with other datasets 
with a better spectral resolution, but not in a number sufficient 
to fit the calibration relation. For this reason we choose to first 
calibrate the archive with this dataset. It was entered in the Hy- 
percat Fits Archive (HFA) with the identification L1985BURN. 

The first step was to assess the photometric quality of these 
templates. The B and V magnitudes were integrated on the 
spectra and compared with the measurements converted from 
the Tycho-2 catalogue. We found a RMS difference on B-V of 
0.05. The difference with the Tycho-2 color is used for affect- 
ing a weight to each of these template observations. 

The determination of the calibration relation is done af- 
ter the instrumental spectra from our archive are convolved to 
the resolution of L1985BURN. The original reference does not 
precisely document this resolution which varies with the origin 
of the observations, and we had to determine it. A first guess 
was obtained by fitting Gaussians on strong Balmer lines and 
this was fine-tuned afterwards to minimize the residuals on the 
strong lines. It was necessary to adopt slightly different resolu- 
tions for the blue and red parts of the spectra. 

The model adopted for the primary calibration relation is: 

log{F phys (\)) = log{F inst (X)) +log(F( nst (X, run)) 
+a(\)Airm + b(X)Airm 2 
+c(\)E B _ v + d(X)E' B _ v (1) 

Fphys is the physical flux reduced above the atmosphere. 
The two right-hand terms \og(F inst (A) and log(F/ n;st (A, run)) 
represent the instrumental response. The first is stable over all 
the observing runs while the second contains the variable part. 
We separated the two terms because the second could not be 
determined for all the observing runs: when the number of cal- 
ibrated observations was less than 10 the second term was not 
computed and assumed to be 0. The run-dependent instrumen- 
tal response was determined for 9 runs out of 21, and a sec- 
ondary calibration (next sub-section) was also obtained. 

The next two terms depend on the airmass (Airm). The first 
of them accounts for the clear atmosphere absorption and mean 
aperture effect due to the atmospheric extinction. The second 
partly corrects the "curvature" of the spectrum due to the re- 
fraction (the length of the atmospheric spectra is proportional 
to tan z, z being the zenithal distance). 

The two last terms are parameterized by a B-V color ex- 
cess, defined as the difference between B-V determined on our 
spectra and the Tycho-2 B-V. 

It is not possible to measure directly the B magnitude on 
our spectra since the standard B band extend farther blue-ward 
than their limit. Hence it was necessary to "calibrate" a differ- 
ent color, by fitting a color equation as it is usually done for 
converting between photometric systems. This procedure has 
the disadvantage that it formally depends on the actual spec- 
tral energy distribution of the star and hence may introduce a 
bias linked to the atmospheric parameters of the stars. How- 
ever, such effect may be estimated from the residuals of the 



5 

fit of the color equation, and for small color terms (ie. if our 
internal B is close to the Johnson B) it can be neglected. 

Since this conversion was necessary, we used this opportu- 
nity to define 2 B bands (b\ and b 2 , converted to Johnson scale 
as Bi and B 2 ) giving two estimates of the color excess. The 
difference between these two estimates gives an approximate 
measure of the residual curvature of the spectra due to the at- 
mospheric dispersion. 

E(B - V) = (B - V) T ycho - (B 1 - V) Elodie (2) 

E(B - V)' = B 2 -B x (3) 

The B\ and B 2 bands were respectively defined as the in- 
tegral of the flux distribution between AA = 450 — 470 nm and 
AA = 437 - 450 nm. 

The color equation was determined with the 2238 spectra 
of L1985BURN, using ordinary least squares and iterative re- 
jection of outliers, as: 

h = 14.830 + l.OOLBi - 0.337(5 - V) (4) 
b 2 = 15.195 + 1.003S 2 - 0.081(B- V) (5) 

The RMS residuals from these fits are about 0.025 mag. 
The magnitudes and colors are practically measured on the in- 
strumental spectra (not physical flux), and hence the measured 
excess is biased due to the instrumental response and depends 
on the airmass. Though it does not change the fit to Eq. 1, in or- 
der to have independent functions we corrected the excess from 
the effect of instrumental response and atmospheric extinction 
using an a posteriori fit between corrected and un-corrected 
excess. 

The flux calibration relation was fitted on 369 spectra of 
stars in common between our archive and L1985BURN. The 
different functions of the wavelength were finally smoothed by 
fitting polynomials in order to erase the residual effects of res- 
olution miss-matching which results in spurious irregularities 
in these functions (in particular around H 7 and G-band). 

4.2. Secondary calibration 

The primary calibration allowed us to determine the run- 
dependent response for only half the runs (but more than half 
the spectra since the number of stars is not equal in each run). 
By inter-comparing spectra of the same stars for which several 
observations were available it is possible to make a secondary 
calibration, using the runs where the run-dependent response 
was available as "template". 

This exercise allowed us to detect a significant change in 
the response of the instrument between the observations prior 
to 1995 and those done after. This change seems to be due to 
a modification of the differential refraction corrector at the end 
of year 1994. In addition, for two other runs the number of 
inter-comparisons was large enough to adopt a reliable run- 
dependent corrections. The run-dependent correction account 
for a RMS photometric effect of about 0.025 mag. For 9 runs 
no such correction is available at present. 
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4.3. Correction of the instrumental response 

Up to this point the flux calibration can be assumed to prop- 
erly correct the effect of atmospheric absorption and to cor- 
rect at the best the aperture effect due to the differential refrac- 
tion. However, the instrumental response was evaluated using 
low-resolution spectra and had furthermore to be smoothed. 
At variance with the other functions which presumably vary 
smoothly with wavelength, the instrumental response may be 
more chaotic. 

Therefore, the last step of the calibration process consists 
in using higher spectral resolution (R w 1000) to correct the 
instrumental response. These reference spectra are taken from 
the following dataset, referenced by their HFA identification: 

- L961 1 1KP1 and L961 1 1KP2. The ELODIE database con- 
tains 191 spectra of 144 stars from the Jones library (Lei- 
therer et al. 1998) (R = 2800). This library cover two nar- 
row wavelength ranges, one in the blue part below A = 
451 nm the other between A A = 478 — 547 nm. This is 
the highest resolution reference for stellar spectra yet avail- 
able for comparison. This library is not corrected for the 
aperture effect due to differential refraction and is not of 
photometric quality. After re-calibrating the "color" of this 
spectra (defining special color-band as we did for our spec- 
tra for the primary calibration), the mean agreement with 
our spectra is satisfactory. 

- L1996SERO. The ELODIE database has 8 spectra of 6 stars 



in common with Serote Ross et al. (1996). These spectra are 



from CFH (R=650) and OHP (R=4000) observations. They 
intersect the red half of our spectral range. The comparison 
with the L1996SERO spectra, after the spectral resolutions 
are matched by convolving our spectra with a Gaussian, 
show pixel-to-pixel variations of the order of 2% (rms) and 
low-frequency variations of typically 5%. These figures are 
consis tent w ith the comparison perfor med b y Serote Ross 
et al. ( |1996| ) with the Silva & Cornell ( |1992| ) library. 

- L1984JACO. We have 2 spectra for 2 stars (HD094028 & 
HD028099) in common with the Jacoby library (R=1200) 
which cover the full wavelength range. The two spectra 
compared with L1984JACO reveal a significant gradient 
(the difference over the whole range of wavelength are re- 
spectively 29% and 14% for the first and second spectrum) 
and a pixel-to-pixel variation of 0.5% rms. The B-V color 
integrated from the L1984JACO spectra disagrees with the 
Tycho-2 value in amounts that explain the differences with 
our spectra. 

- L1987KIEH. We have 13 spectra of 7 stars from the dataset 
of Kiehling (1987) (R=800). 

- L1994DANK. We have 17 spectra from 14 stars in common 



with Danks & Dennefeld (1994) (R=1500). The red half of 



our wavelength range intersects with L1994DANK. 



- L1983GUNN. The Gunn & Stryker (1983) library (R=500) 



counts 5 stars in common with our archive. 
- L1984NSOA. Finally we compared our 7 solar spectra with 
the Solar flux atlas (Kurucz et al. 
stored in HFA with R=10000). 



1984) (R=500000, also 



We made all the individual comparisons between our spec- 
tra, convolved by a Gaussian to match the resolution, and these 
references. Because the broad-band variations are probably due 
either to random photometric errors in our archive (certainly 
not systematic), correction of the interstellar extinction (it was 
not always possible to de-correct the templates) or errors in the 
templates, we subtracted a low-degree polynomial (1 to 3 de- 
pending for which dataset) to these comparisons. 

The individual comparisons (see Fig. |3|) were then com- 
bined and produced a mean residual convolved to a resolution 
of 2 nm FWHM (each order of the original spectrum covers 
about 4 nm and no correction of the instrumental response at a 
lower scale is desired). This was used to correct the instrumen- 
tal response. 

We also made an un succes sful attempt to compare our 
archive with the Pickles ( |1998| ) library (L1998PICK): all the 
spectra were compared to the nearest (in the T c ff - log g - 
[Fe/H] space) template in L1998PICK. But the mapping of 
the parameters space in L1998PICK is not dense enough and 
the residuals are too affected by the distance in the parameter 
space to permit a useful analysis. We did not attempt to inter- 
polate the spectra within the Pickles library. 

5. Quality control 

To check the reliability of our data and reduction we used 
the following tests: (1) analyze the inter-comparison of re- 
peated observations and analyze internal consistency, (2) com- 
pare fully reduced spectra with other libraries and (3) compare 
Lick indices for stars with published measurements. 

There are a priori three possible sources of discrepancy to 
search: (1) errors on the instrumental response and/or correc- 
tion of atmospheric extinction are likely to result in low and 
intermediate spatial frequencies variations. (2) inaccurate sub- 
traction of the diffuse light may be diagnosed as significant 
mismatch of the strongest lines and (3) inaccurate correction 
for the blaze will lead to oscillations with a one-order period. 



5.7. Internal consistency 

5.1.1. Multiple observations of the same star 

There are in total 358 repeated observations for 104 stars in 
the archive. After giving a double weight to the observations 
in the runs for which the run-dependent response was com- 
puted from primary calibration and half weight to the obser- 
vations done with an airmass larger than 1 .5 (their photomet- 
ric quality is lower because of the differential refraction), we 
find a global RMS dispersion for the pairwise comparison of 
about 3.5%, see Fig. [|. Without weighting, the dispersion is 4% 
and the difference was a hint for the lowest photometric qual- 
ity of high airmass observations, hence leading to suspicion 
about the refraction corrector. Since it is a pairwise compari- 
son, the mean photometric error on individual observations is: 
3.5 /-s/2 = 2.5%. (Note that this error does not account for pos- 
sible systematic errors affecting uniformly the whole archive, 
it will be discussed in Sect. 5.2 & 5.3). 
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Fig. 3. The mean photometric residuals of the 9 libraries used for the correction of the instrumental response (before this correc- 
tions). The ordinate is the mean difference (Elodie - reference library) in units of relative physical flux (normalized to 1 at 555 
nm). Each dataset is shown on a different line, each being arbitrary y-shifted by 0.05. From bottom to top they are: Jones (blue 
and red), Serote Roos et al. (OHP and CFH), Jacoby, Kiehling, Danks & Dennefeld, Gunn & Stryker and the Sun. For the purpose 
of the display these spectra have been slightly smoothed. The top line is the combination of all the individual comparisons which 
were used to compute the correction to the instrumental response. 
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Fig. 4. Photometric precision deduced from the comparison of multiple observations of the same star. The ordinate are the rms 
of the 358 pair comparisons (ln(Flux 2 /Fluxi)) (the spiky appearance is due to the logarithmic comparison). The higher curve 
the the total RMS while the lower one is after third degree polynomials were subtracted from each individual comparison: This 
disentangles the effect of photon noise and uncertainty on the long range flux calibration (due to the differential refraction). To 
obtain the order of magnitude of the photometric error on single observation, the ordinates should be divided by \/2. 



Most of the photometric errors lie in residual "curvature" 
of the spectra. After subtracting a third degree polynomial to 
each pairwyse comparison, the RMS dispersion becomes 0.7%, 
corresponding to a mean photometric precision of 0.5%, in full 
agreement with the mean S/N ratio. No hint for a modulation 
at the scale of one order of the echelle spectra is found. 

5.1.2. Modeled spectra 

Random errors on the strengthening of the spectra are probably 
in average smaller than 0.5% because otherwise they would in- 
crease the dispersion of the pairwyse internal comparison (after 
the third degree polynomial is subtracted). 



The new version of the TGMET program that will be de- 
scribed in a separate paper allowed us to generate interpolated 
spectra for a given position in the (T ff , log g , [Fe/H]) space. 
The difference between the spectra in the archive and these 
modeled spectra can potentially uncover possible random er- 
rors on the strengthening. 

Therefore we computed and averaged the differences to the 
modeled spectra for all the spectra in the archive. The mean 
residuals were searched for a modulation at the scale of one 
order, either in the spectrum or in its Fourier transform. We 
did not detect any signal. We made simulations that told that a 
0.5% (RMS) modulation would probably have been detected, 
but a 0.2% modulation would not. 



8 



Ph. Prugniel & C. Soubiran: Database of stellar spectra 



5.2. External comparisons 

The comparisons between our physically calibrated spectra and 
other libraries were used in the previous section to compute 
the correction to the instrumental response. The goal was to 
determine the variation of the response at the scale of 1 to 5 nm 
that was not constrained by the primary calibration. Therefore, 
by construction, the mean external comparison is fiat (ie. may 
only varies smoothly with the wavelength). A careful analysis 
of the individual comparison could assess the level of stability 
of the instrumental response. 

The patterns in the correction to the instrumental response 
were consistently found in the different individual comparisons 
and no significant effect was found after this correction was 
applied. However, the scarcity of the available comparisons did 
not allow us to measure a photometric precision. 

5.3. Comparison of Lick indices 

The Lick indices (see eg. W097) measure equivalent widths of 
features defined in bands of 2 to 5 nm (hence of the width of 
about 1/2 or 1 order of the echelle spectra). Their measurements 
are then sensitive to both incorrect strengthening of the spectra 
and errors in the instrumental response. 

The spectral range of the Elodie archive allowed us to mea- 
sure most of the Lick indices. Only the blue CN indices are 
missed as well as the H<5 indices (H5a and HSf) defined in 
W097. 



Table 2. Comparison of Lick indices: Col. 1: Name of the 
Lick index Col. 2: Number of spectra used for the compari- 
son Col. 3: Mean value of the index Col. 4: Mean shift (Elodie) 
- (W097) Col. 5: RMS dispersion of the comparison Col. 6: 
Majorant of photometric error 



We used the definitions from W097 or Cardiel et al. fll998| ) 
to measure the indices for all the spectra in the archive. The 
results are reported in Table 1 given in the electronic version 
only. 

The measurement errors computed from the photon noise 
(see eg. Prugniel et al 2001 1 is in the present case negligible 
because of the extremely high S/N per A. The main source of 
error is due to the photometric calibration and may be system- 
atic. The comparison the Lick indices measured in W097 con- 
strain the magnitude of these systematic errors. 

The archive has 102 stars which have been previously mea- 
sured in W097. The comparisons are reported in Table 2. 

We do not find systematic effects that would result from 
errors in the instrumental errors and the residual dispersions 
are compatible with the errors in the W097 measurements (the 
internal error due to the photon noise in the archive is totally 
negligible). 

The systematic shifts are always smaller than the dispersion 
except for Ca4455 (where it is equal) and for Fe4383 the shift 
may be significant. However, we cannot attribute this to an er- 
ror in our measurements and in particular we note that the shift 
for Fe4383 is similar to the one found by W097 when they 
compared the original Lick measurements with the spectra of 
the Jones library (note that our flux calibration is not fully inde- 
pendent from the Jones library since this latter was used among 
others to correct the instrumental response). The Ca4455 index 
is not in the range of the Jones spectra and we cannot comment 
of the shift we found. Fig. § presents the distributions of the 
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differences (us - W097) for two indices: Mg2 and H-ja- The 
first one is very important for extragalactic studies, it is defined 
on a broad region (53 nm between the extrema of the red and 
blue side-bands). Though it is potentially sensitive to the uncer- 
tainties on the flux calibration described above, the agreement 
between us and W097 is excellent. The second index presented 
in Fig. 5, Hja, is narrower (14 nm) but it is located in the blue 
region of our spectra were our calibration is probably less ac- 
curate. There is a systematic shift with respect to W097, but 
the rms residual is fully accounted for by the errors on W097. 

6. Data format and public access 

The spectra are stored in HFA in the form of FITS files con- 



taining 4 extensions (Gr0sbol et al 1988): 

The first extension is a ID image header data unit (HDU) 
containing the spectra calibrated in wavelength expressed in the 
rest-frame of each star and sampled between 410 and 680 nm. 
The high resolution version is normalized with a step between 
the sampling points of 0.005 nm and the low resolution version 
is given in physical flux with a step of 0.020 nm. 

Another image HDU (EXTNAME=NOISE), with the same 
sampling, contains the noise attached to the values in the pri- 
mary HDU (in the unit of this spectra). Pixels are independent. 

The other extensions are the flux calibration relations which 
enable to convert the instrumental flux into physical or normal- 
ized flux. They are binary table extensions and the relations 
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Fig. 5. Comparison of Lick indices Mg2 and Hja with W097. 
The abcissae are the difference (Elodie - W097) between the 
Lick indices. The top panels show the distribution as a function 
of the value in W097, and the bottom panels the histograms of 
the differences. 

are coded as spline representations. The details are given in the 
Pleinpot user's manual^], and the conversions can be done on- 
the-fiy with the pipeline attached to HFA. 

The physical flux calibration (EXTNAME=FCAPHY) and 
the normalized flux relation (EXTNAME=FCANOR) have a 
spectral resolution of about 2 nm. 

Each header, which includes the keywords created dur- 
ing the observations and TACOS reduction, have been 
completed with useful information on the star : spec- 
tral type, Tycho-2 V and B-V, absolute magnitude My and 
(T c ff , log g , [Fe/H]) with their reliability flag, and v r . The 
nominal FWHM of the spectrum can also be retrieved in the 
header. 

In addition to the spectra of the stars, the archive contains 
the different intermediate files used along the data processing. 
They are described in the electronic version of the paper. 

7. Conclusion 

We have presented an archive of 908 spectra of 709 stars at 
two resolutions, R= 10000 and R=42000, and in the wavelength 
range AA = 410 — 680 nm. The high resolution version is 
archived normalized to its pseudo-continuum while the low res- 
olution is given in physical flux (above the earth atmosphere) 
normalized to one at A = 555 nm. 

The mean S/N per pixel in the high resolution version of the 
archive is 150 (S/N w 700 per 0.1 nm), the photometric preci- 



sion of the strengthening of the orders of the echelle spectra is 
better than 0.5 %, the precision in bands of about 10-20 nm is 
0.5 % and the overall (broad-band) precision is 2.5 %, limited 
by the effect of differential refraction. 

The stars in this archive cover a large range of the parame- 
ter's space, though the mapping is still crude for hot and very 
cool stars. 



3 http://www-obs.univ-lyonl.fr/~prugniel/cgi- 
bin/pleinpot/pleinpot.html 



This archive supersedes the TGMET ( 1996 ) library and is 
being used for automatic determination of the stellar parame- 
ters (T e ff , log g , [Fe/H]) by direct comparison of spectra with 
interpolations in the archive. The new version of the TGMET 
program is still under development and will enableto determine 
the parameters from any R> 10000 spectra, independently of 
the detailed observational setup. It will be presented in a future 
paper. 

The archive will also be used to synthesize extragalactic 
stellar populations and to use these high resolution synthetic 
spectra as kinematical templates to determine simultaneously 
the parameters of the stellar population and the line-of-sight 
velocity distribution. 

The archive itself will be completed with new observations 
and spectra from the OHP archive with the goal of improving 
the mapping of the parameters space and the quality of the pho- 
tometric calibration. 
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